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We experimentally and numerically demonstrate deterministic creation and manipulation of a pair 
of oppositely charged singly quantized vortices in a highly oblate Bose-Einstein condensate (BEG). 
Two identical blue-detuned, focused Gaussian laser beams that pierce the BEG serve as repulsive 
obstacles for the superfluid atomic gas; by controlling the positions of the beams within the plane of 
the BEG, superfluid flow is deterministically established around each beam such that two vortices 
of opposite circulation are generated by the motion of the beams, with each vortex pinned to the in 
situ position of a laser beam. We study the vortex creation process, and show that the vortices can 
be moved about within the BEG by translating the positions of the laser beams. This technique 
can serve as a building block in future experimental techniques to create, on-demand, deterministic 
arrangements of few or many vortices within a BEG for precise studies of vortex dynamics and 
vortex interactions. 


I. INTRODUCTION 

Quantized vortices are central features of macroscopic 
quantum coherent phenomena urn, and serve as robust 
indicators of the dynamics of quantum fluids. Through 
the development of progressively better experimental 
techniques for controlling and observing the motion of 
vortices in these systems, new possibilities open up to 
better understand the physics of vortices in a quan¬ 
tum fluid. In turn, new techniques for controlling the 
states and properties of quantum fluids may soon fol¬ 
low. For example, in type-II superconductors immersed 
in a strong magnetic field, the motion of magnetic flux 
vortices is largely responsible for energy dissipation and 
limiting the temperatures at which superconductivity in 
these materials occurs [3] ; by pinning or otherwise arrest¬ 
ing the motion of vortices, superconductivity may persist 
at higher temperatures and magnetic fields than if the 
vortices were free to move. New research on the param¬ 
eters that influence vortex motion may thus help extend 
the temperature and magnetic field ranges of such su¬ 
perconductors [5]. In superfluid helium, vortex-vortex 
interactions and Kelvin waves on vortex lines also dissi¬ 
pate energy; new methods to measure and better under¬ 
stand the dynamics of vortices in a superfluid might shed 
new information on quantum turbulence energy dissipa¬ 
tion mechanisms in these systems [6]. Similarly, in two- 
dimensional (2D) quantum turbulence in highly oblate 
Bose-Einstein condensates (BECs) [7], vortex motion is 
a key indicator of the hydrodynamic state of the BEG; 
studies of vortices and their motion in these systems may 
aid in the development of a deeper understanding of 2D 
quantum turbulence @19]. To better understand these 
and other quantum fluid phenomena, it is essential to 
construct a detailed understanding of vortices, their dy¬ 
namics, and their interactions under a wide variety of 
conditions and in different systems. 

Dilute-gas BECs are well suited to experimental, the¬ 


oretical, and numerical studies of vortices, and to direct 
quantitative comparison between experimental and nu¬ 
merical results that enable rapid development of our un¬ 
derstanding of vortices and superfluidity in BECs. With 
the publication of over 100 articles to date covering ex¬ 
perimental work with vortices since 1999, the experimen¬ 
tal techniques available for creating and studying vortices 
in BECs are by now numerous and widely varied mm- 
In most of these studies, the specific positions of each 
vortex within the BEG are not of primary concern, and 
vortex creation is not deterministic at the single-vortex 
level. Nevertheless, deterministic vortex creation is pos¬ 
sible, and some experiments have indeed focused on de¬ 
terministic creation mechanisms and even pinning of a 
singly quantized or multiply quantized vortex or estab¬ 
lishment of a persistent current within a BEG p!2H24] . 
Vortices in BECs are also large enough to be detected 
optically, either with in situ techniques involving multi- 
component BECs m or single-component BECs |25], or 
most commonly with a now-standard method of imaging 
the BEG after a short period of ballistic expansion [26] . 
Various techniques enable the dynamics of vortices to be 
measured, either in real-time 1271 HE] or after determin¬ 
istic vortex-dipole nucleation [29] for single-component 
BECs, or with in situ imaging for multi-component BECs 
[12]. Vortex dynamics and pinning are also strongly af¬ 
fected by the presence of time-dependent potential bar¬ 
riers [30] , and even arrays of vortices can be pinned in a 
rotating frame by laser beams [31]. 

However, despite such a proliferation of experiments 
with vortices in BECs, there are still no published ex¬ 
perimental techniques that have demonstrated the flex¬ 
ibility of constructing arbitrary vortex configurations in 
a BEG, such that the location, vorticity, and circulation 
direction for each individual vortex can be presicribed de¬ 
terministically. Such a technique would be an important 
component of precision studies of vortex dynamics and 
interactions, and in manipulating and studying states of 
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a quantum fluid. Especially under conditions where ex¬ 
isting single-shot imaging techniques are used, the abil¬ 
ity to specify vortex distributions on-demand would open 
up new possibilities for studying vortex dynamics and in¬ 
teractions. To address the challenge of developing such 
a technique, we present here a method of deterministi¬ 
cally generating and manipulating two singly quantized 
vortices of opposite circulation in a highly oblate BEC, 
showing results of an experimental study and numerical 
simulations. We expect that this technique is scalable 
well beyond two vortices, and in principle may be used to 
create designer distributions of many vortices; this possi¬ 
bility is numerically explored in a companion article [32] . 

Below, we first present the conceptual foundation for 
our technique, followed by a description of our experi¬ 
mental setup and the beam parameters used in our study. 
To illustrate the vortex creation process described above, 
we then show results of numerical simulations of the 2D 
Gross-Pitaevskii equation for conditions closely approxi¬ 
mating those of our experiment. Einally, we discuss var¬ 
ious aspects of the manipulation of vortices within the 
BEC. 


II. CONCEPT 

Our technique generally proceeds as follows. We begin 
with two blue-detuned laser beams that pierce a station¬ 
ary BEC at positions that can be controlled electroni¬ 
cally. The beams are initially co-located, then simulta¬ 
neously and at a constant speed move together across the 
BEC, say in the x-direction, while also moving slightly 
apart in the ^-direction before coming to rest at posi¬ 
tions for which the beams are fully separated. The laser 
beams serve as repulsive obstacles for the atomic super¬ 
fluid, pushing fluid out of their way as fluid simultane¬ 
ously fills the space vacated by the beams during their 
motion. The beam motion results in macroscopic flow 
within the BEC, even for conditions where the speed of 
the laser beams is far too slow for vortex dipole nucle- 
ation and shedding into the bulk superfluid. 

As the beams move through the BEC, they eventually 
separate and become two moving obstacles, allowing the 
flow that was occupying the space behind the combined 
beams to merge with the portion of the fluid that was 
being pushed out of the way ahead of the beams. Eor a 
suitable range of beam trajectory parameters, the now- 
continuous fluid flow around each separate beam corre¬ 
sponds to a net winding of the BEC’s quantum phase 
of ±27r. In other words, at the moment of beam separa¬ 
tion, two opposite-circulation regions of persistent super¬ 
fluid flow are trapped around each laser beam. Equiva¬ 
lently, this process can be seen as the simultaneous nu- 
cleation and pinning of two oppositely charged vortices, 
with one vortex pinned to each beam. Subsequent motion 
of the laser beams can transport vortices to other loca¬ 


tions within the BEC, reminiscent of the use of lasers as 
optical tweezers for material particles, and allowing for 
deterministic vortex creation and location. To emphasize 

(i) the need to use laser beam pairs in our new method, 

(ii) the two different linear trajectories that the beams 
take when moving through the BEC (described in detail 
below), and (hi) the vortex pinning and manipulation 
capabilities provided by the beams, we henceforth refer 
to this dual-vortex creation and manipulation method as 
the “chopsticks” technique. 

III. EXPERIMENT 

We create BECs of ^^Rb in the 5^5'i/2 |E = 1, mi? = 
— 1) state in a highly oblate harmonic trap. The trapping 
potential is a combination of a magnetic time-averaged 
orbiting potential (TOP) magnetic trap [33] with an 
axis of symmetry in the vertical (z) direction, and a 
red-detuned 1090-nm laser beam that propagates hor¬ 
izontally along the x direction and is tightly focused 
and provides tight confinement in the z direction. The 
harmonic trapping frequencies of this hybrid trap are 
{ujrl^TT^ujzl^Tr) = (8, 90) Hz in the radial and vertical di¬ 
rections, respectively. With no laser beams present other 
than the 1090-nm laser beam used for the hybrid trap, 
our BECs have Thomas-Eermi radii of Rr ^ 50 pm in 
the radial (x^y) plane, corresponding to condensates of 
about Nc ^ 2 X 10^ atoms and a chemical potential of 
Po ^ 8 Eurther details of our BEC creation methods 
can be found in Refs. [71|29|. 

The primary tool used in the chopsticks technique is 
a pair of blue-detuned focused laser beams propagating 
along the z direction. The beams serve as repulsive bar¬ 
riers that penetrate the BEC and whose positions can be 
translated within the plane of the BEC. Eigure 1 shows 
a schematic diagram of the optical system that is used 
to control the motion of these two beams; a detailed de¬ 
scription of the optical system is provided in the figure 
caption. The chopsticks beams are initially co-located 
at position {xi, yi) = (—10, 0) pm relative to the cen¬ 
ter of the trap. Prior to the final stage of evaporative 
cooling that creates a BEC, each beam is turned on to 
full power, corresponding to a peak repulsive potential 
of ^ 0.8 po- By creating the BEC after the beams are 
in place, we minimize excitations of the BEC that would 
otherwise result from turning on the beams after creating 
a BEC, and also damp out BEC sloshing in the harmonic 
trap that can otherwise decrease the effectiveness and re¬ 
peatability of vortex creation. Eigure 2(a) shows an in 
situ image of the BEC with the two chopsticks beams 
turned on at full power at their initial locations. 

To generate two vortices of opposite charge pinned 
to the chopsticks beams, we simultaneously move both 
beams along linear trajectories at a constant speed for a 
sweep time of C = 0.7 s. Beam Bl, in the upper half- 
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(a) Optical fiber 



FIG. 1. (Color online.) (a) Diagram of the optical system that 
is used to translate two blue-detuned laser beams through the 
BEG. A 660-nm laser beam is brought into the optical sys¬ 
tem with a polarization-preserving single-mode optical hber. 
The beam is collimated by lens LI and then split into two 
paths by a polarizing beamsplitter cube (PBSC); the optical 
powers of the two beams are balanced using a half-wave plate 
(HWP) placed before the first PBSC. Each beam reflects off 
of two mirrors (Mia and Mlb, or M2a and M2b). Each mir¬ 
ror has a piezo-electric transducer (PZT) stack attached to an 
adjustment screw; the two PZTs attached to the mirrors of 
a given beam path control the tilts of the two mirrors about 
orthogonal axes, enabling the two laser beams to be remotely 
steered. After passing through additional HWPs, allowing for 
further optical power adjustment, the two beams are recom¬ 
bined with another PBSC. Lenses L2 and L3 comprise a 6x 
minihcation system in order to increase the angular deflection 
of the chopsticks beams, prior to deflection of the beams by 
mirror M3 and focusing the beams at the BEC by lens L4. In 
the experiments reported here, the beams have 1/e^ radii at 
the BEC of (Jc ~ 19 /xm. During the course of an experimen¬ 
tal run, the two laser beams move along paths PI and P2 in 
the (x, y) plane, as represented in the inset diagram, (b) An 
illustration of the BEC (large gray disk) lying in the {x,y) 
plane. The laser beams are initially co-located and pierce the 
BEC (left dark-gray disk) and move along paths PI and P2, 
which subtend an angle a ~ 66°. The beams end up in dif¬ 
ferent locations (upper and lower right dark-gray disks). The 
size of the laser beams, relative to the BEC Thomas-Eermi 
radius, and the initial and final locations of the laser beams, 
are illustrated to scale. 


plane, moves to a final position of (x/, yf) = (13, 15) /xm. 
Beam B2, in the lower half-plane, moves to a final po¬ 
sition (xj, yf) = (13, —15)/xm. An in situ image of 
the BEC after the beams have been moved is given in 
Fig. 2(b). The beams travel distances of ^ 27/xm at a 
speed of ^ 34 /xm/s, much less than the maximum speed 
of sound for our BECs of cq ^ 1700 /xm/s, and still much 
less than ^ cq/IO, the speed required to nucleate and 
shed a vortex dipole for comparable beam and BEC pa¬ 
rameters m- Once the beams reach their final positions. 


(a) I (c) 





EIG. 2. (a) A 144-/xm-square in situ image of the BEC 

taken along the vertical [z) axis, showing the prohle of the 
~100-/xm-diameter BEC in the {x^y) image plane. The hole 
to the left of center is created by the two co-located chop¬ 
sticks beams turned on at full power prior to moving them 
across the BEC. (b) Same as (a), but with the chopsticks 
beams at their final locations at the ends of paths PI and 
P2. The intensity minima due to the chopsticks beams are 
~ 31 /xm apart in this image, slightly larger than the sepa¬ 
ration obtained when measuring the beam positions directly, 
(c) A 205-/xm-square absorption image of a BEC at the end 
of the chopsticks procedure, after the chopsticks beams were 
ramped off at their final positions and the BEC was released 
from the trap and allowed to ballistically expand a factor of 
~ 1.78 in the radial direction. Vortex cores appear as the two 
small dark holes in the density distribution, and are ~ 65 /xm 
apart in this image, compared with an expected separation of 
~ 55 /xm given the hnal positions of the chopsticks beams and 
the expansion factor of the BEC. 

their positions and powers are held constant for times up 
to 2.5 s, after which their optical powers are linearly de¬ 
creased to zero in 0.25 s. This rate at which the beams 
are linearly ramped off was determined experimentally, 
such that it is fast enough to ensure that the vortices 
would not move significantly far away from the beam po¬ 
sitions as the beam power is decreased, but slow enough 
such that the BEC is not significantly excited by a rapid 
change in the potential well profile. The BEC is then re¬ 
leased from the trap and allowed to ballistically expand 
in order to look for the presence and positions of vortices 
using standard absorption imaging techniques. 

Fig. 2(c) shows a representative image in which two 
vortices were observed after applying the chopsticks pro¬ 
cedure, with each vortex being located near the final posi¬ 
tion of one of the laser beams. In one test of the repeata¬ 
bility of the chopsticks procedure, we created 30 BECs 
and subjected each to the procedure described above. 
In 26 (86%) of the cases we found two vortices, with 
each located at a position that closely matches the fi¬ 
nal beam position, relative to the BEC center. In the 
remaining cases, either one vortex was observed in the 
BEC, or an extra unpinned vortex was observed in ad- 
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dition to the two vortices that were located at positions 
matching those of the chopsticks beams. A likely expla¬ 
nation for these cases is the presence of spontaneously 
formed vortices [34]. Such vortices can alter the back¬ 
ground flow profile of the BEC such that the necessary 
conditions for vortex nucleation and pinning parameters 
are not achieved, or might annihilate with one of the vor¬ 
tices created by the chopsticks technique, or could remain 
present throughout the chopsticks mechanism of vortex 
creation and be visible as a third vortex. 

IV. NUMERICAL SIMULATIONS 

To better understand the mechanisms at work in the 
chopsticks technique, we turn to numerical simulations 
of the Gross-Pitaevskii equation (GPE) [35]. We first 
assume a three-dimensional wavefunction of the form 

where the dynamics of the 2D macroscopic wavefunction 
' 0 (x, t) are assumed to closely approximate the dynam¬ 
ics of our actual three-dimensional (3D) system. This ap¬ 
proximation is appropriate for our highly oblate system, 
particularly because for our conditions, vortices remain 
aligned along the ^ direction without undergoing signifi¬ 
cant tilting or bending [36| . The constant Iz is an effective 
length scale for the axial thickness of our BEG, described 
below, and is not necessarily equal to the harmonic os¬ 
cillator length for the z direction. 

With this 2D approximation, the 3D GPE with a di¬ 
mensionless phenomenological damping constant 7 can 
be reduced to 

{i - + V +52d|V’|V, 

where Wt = is the potential energy of a 

2D harmonic trap of frequency and m is the mass of 
^^Rb. The potential due to the chopsticks beams is given 
by 

2 

Vc{x,y,t) = Uo'^exp{-2[{xi,yi)j + {v^,, Vy)jt]'^/a^} 
i=i 

where Uq = 0 . 8 /io is the maximum repulsive potential en¬ 
ergy of each beam during the sweep, {xi^yi)j is the initial 
position of beam j, {vx, Vy)j is the sweep velocity vector 
for beam j, and cTc = 19 /im is the radius of each chop¬ 
sticks beam. In order to mimic the effects of non-zero 
temperature and the damping provided by a small ther¬ 
mal cloud of atoms, 7 is set to 0.003 for all simulations 
reported here mim- 

We use a 2 D interaction parameter ^ 2 D = , 

where a is the atomic 5-wave scattering length for our 
atoms. Eor the 2D GPE simulations, we use an effective 



0 , 1 0 ^ 2jt 0 , 0.36 

n/n^ax ^ v/Co 


FIG. 3. (Color online.) Results of 2D GPE numerical simu¬ 
lations of the chopsticks technique for simulation parameters 
described in the text. Each row represents the state of the 
system at a time t during the sweep, with the sweep time t/ts 
indicated in the upper left corner of the images of the left 
column. The final row represents the state of the system at 
a hold time of 0.5 s after the end of the sweep, during which 
time the chopsticks beams’ positions and powers remain un¬ 
changed from their final values. Columns one through three 
show (respectively) for each time step 120-/im-square regions 
of the density profile n relative to the maximum number den¬ 
sity Umax, the quantum phase profile 0, and the velocity pro¬ 
file relative to the maximum speed of sound co = 
for the ground state BEC without chopsticks beams present. 
Column four shows a magnified view of the velocity profile, a 
61-/im-square region centered on the point midway between 
the chopsticks beams. Arrows indicate the local velocity of 
fluid flow at the tail of each arrow. For clarity, in all columns, 
points for which n/umax < 0.0005 are assigned the color black. 
The colored circles in the second column are described in the 
caption of Fig. 4. 
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length parameter = 0.53/rm, an effective radial trap 
frequency of = 0.84 where Ci;^/27r = 8 Hz as in 
the experiment, and an effective number of atoms equal 
to 3.2 X 10^. These numbers are chosen such that when 
the BEC is in its ground state without chopsticks beams 
present, the density at the center of the BEC in the 2D 
approximation matches that of a 3D Thomas-Eermi ap¬ 
proximation for the actual experimental parameters, and 
the 2D Thomas-Eermi radius of '0(x,^) is also the same 
as that of the 3D Thomas-Eermi approximation. This ap¬ 
proach enables us to compare the beam velocities of the 
simulation to a speed of sound that is the same as calcu¬ 
lated for our 3D experimental parameters in the Thomas- 
Eermi limit, Co = permits direct comparison 

of beam radii and positions in the simulations with the 
equivalent parameters of the experiment. 

Using a split-step routine, we simulate the dynamics of 
'0(x, y, t) on a 120-/im-square area using a 512 x 512 grid, 
giving a grid point spacing of 0.23 /im that is comparable 
to the bulk healing length of ^ 0.28 ym calculated for 
our experimental parameters. After finding the ground 
state of the system in the presence of the static chopsticks 
beams at position {xi,yi)i = {xi^yi )2 = (—10, 0) /im, the 
beams are swept to their final positions (13, ±15) /im in a 
sweep time ts = 0.7 s, as in the experiment. Einally, with 
the beams remaining at their final positions and at their 
full powers, we propagate the simulation for an additional 
hold time of 0.5 s. The results of this simulation are 
shown in Eig. 3. 

We see from the simulation results of Eig. 3 that when 
the two beams start moving, fluid flow is induced in 
the BEC. Eor example, at t/tg = 0.3, the velocity pro¬ 
file shows the highest fluid flow speeds at the sides of 
the moving barrier, corresponding to the regions of the 
largest phase gradients seen in the phase profile. As the 
chopsticks beams eventually separate, between t/tg = 0.7 
and t/tg = 0.8, a channel opens up between the beams, 
allowing the atomic superfluid in the regions in front of 
and behind the beams to merge. In order for to be 
single-valued and continuous in a path around each laser 
beam, the two merging regions of superfluid must adopt a 
continuous phase profile in the channel between the two 
beams. Whether or not the path around a given laser 
beam encloses net circulation will depend upon how the 
regions merge and a phase gradient is established [39] . 
As shown in Eig. 4, prior to beam separation and fluid 
merging, there is a phase difference of approximately 
1.757r established between the regions in the front and 
rear of beam travel; once the the beams separate and the 
front and trailing regions merge, the phase relaxes to an 
azimuthal phase winding of 27r (—27r) around beam B1 
(B2), indicating the presence of a positively (negatively) 
charged vortex pinned by beam B1 (B2). 

The phase winding around each beam at the end of 
the chopsticks process depends on the fluid flow speed 
induced by beam motion. If, immediately prior to merg- 
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FIG. 4. (Color online.) For each time step t/tg represented 
in Fig. 3, the quantum phase (j) around the path denoted by 
the colored circles in the phase profiles of Fig. 3 is shown 
here as a function of the angle 0 around the circle. The point 
^ = 0 is taken to be at the midpoint between the two laser 
beams, and 0 increases in a clockwise direction, as indicated 
by the inset figure showing the phase profile at t/tr = 1.0. 
The plots of 0 vs. ^ in this figure are labeled with the values 
t/tg of the corresponding rows of Fig. 3, and the colors of 
the corresponding rings that indicate the points for which 0 
is shown. At points where the density is negligible (i.e., for 
black-colored values in the phase profiles of Fig. 3), values for 
(j) are not shown. This plot shows the gradual establishment 
of a clockwise 2tt phase winding around the lower laser beam. 
At time t/tg = 0.7 (green line), immediately before the beams 
separate, a phase winding close to 27r has been established. 
At t/tg — 0.8 (pink line), just after the beams separate, a 
continuous loop of fluid is seen to exist around each beam 
(see first column Fig. 3), with a phase winding of around 
each beam. For ease of comparison of the phase at different 
times t/tg^ each plot of 0 vs. ^ shown here has a constant 
value added such that </> is continuous from ^ = 0 to ^ = 27r 
and the minimum value of 0 over this range is zero. 


ing, a phase difference precisely equal to tt has gradually 
built up between these two regions (due to, for instance, 
slower beam velocities than used in our experiment), the 
system will support the temporary existence of a dark 
soliton between the two beams [10]. If the phase differ¬ 
ence between the front and trailing regions is less than tt, 
as is the case in Fig. 5 for tg = lA s, no vortices will be 
created or pinned. If that phase difference is somewhat 
greater than tt and less than 37r, the system will relax to 
one with a positively charged vortex pinned to beam Bl, 
and a negatively charged vortex pinned to beam B2, as 
seen in rows 2 through 4 of Fig. 5. Numerically, the pres¬ 
ence of a small but non-zero damping constant 7 aids 
this relaxation process, but is not essential to the cre¬ 
ation of vortices or their pinning at a laser beam. For 
even larger phase differences, induced by faster sweeps, 
larger amounts of vorticity are created, but the vortices 
are not stably pinned by the laser beams for the con¬ 
ditions studied. The numerical results shown in Fig. 5 
demonstrate the robustness of the chopsticks method to 
substantial variations in the sweep time. 






6 



FIG. 5. (Color online.) Each column shows 61-/im-square 
regions of the phase profile (t){x^y) centered on the midpoint 
between the chopsticks beams for different times t/ts, indi¬ 
cated above each column. Sweep times ts are indicated at the 
left of each row. The display colors are mapped to a range 
of 0 to 27r, as in the second column of Fig. 3. Except for the 
variation of ts , simulation parameters are identical in all cases 
to those described in the text. For 0.6 < ts < 1.2, two vortices 
are generated and left pinned to the beams, an indication of 
the robustness of the chopsticks technique. For ts > 1.2, an 
insufficent phase gradient is established around each beam, 
and vortices are not created. For ts < 0.6, too large of a 
phase gradient is established, and two or more vortices are 
nucleated for each beam. For the conditions studied, these 
additional vortices did not remain pinned to the laser beams, 
as can be seen from the singularities in the phase profiles in 
the last two rows. 


V. MANIPULATING VORTICES 

As an experimental demonstration of deterministic 
manipulation of vortices, after generating and pinning 
two vortices and then holding the beam positions and 
powers fixed for 0.5 s, one of the beams was trans- 



FIG. 6. Experimental demonstration of vortex position ma¬ 
nipulation. After creating and pinning the two vortices, and 
waiting an additional 0.5 s, one of the chopsticks beams was 
returned to its initial position in 1.0 s while the other beam 
remained in place. For the left image and inset, the upper 
beam (Bl) was returned to its initial position, as shown in 
the inset in situ image. The inset for the right image shows 
the result of returning the lower beam (B2) to its initial posi¬ 
tion. In each case, after a beam returned to its initial position, 
both beams were ramped off, releasing the vortices into the 
BEG. The main images of this figure show the corresponding 
positions of the vortices after the BEG was released from its 
trap and imaged following a period of ballistic expansion. 

lated back to its original position over a return time of 
tr = 1.0 s while the other beam remained stationary. 
The beams were then held in these new positions for an 
additional time of 1.0 s. The insets to Fig. 6 show in 
situ images of the BEG with one of the two chopsticks 
beams returned to its initial position after completion of 
the chopsticks procedure. After ramping off the beams in 
0.25 s, and releasing the BEG from the trap, two vortices 
were observed in the BEG at locations that correspond to 
the new final positions of the chopsticks beams, as shown 
in the main images of Fig. 6. 

The chopsticks beams can also be recombined as a 
check that the equal and opposite circulations can be 
cancelled out on-demand. We experimentally examined 
two different methods of returning the beams. First we 
examined sequential beam motion, with one beam re¬ 
turning to the initial position yi) = (—10, 0) jivn in 
1.0 s, followed by a 0.5 s hold during which both beams 
were stationary, followed by the second beam returning 
to the initial position in 1.0 s. We also examined simul¬ 
taneous beam motion, where both beams returned to the 
initial location in 1.0 s. As expected, after returning both 
chopsticks beams to their original position, either sequen¬ 
tially or simultaneously, vortices were repeatedly absent 
from the BEG. From our experimental data, however, it 
is not clear whether vortex-antivortex annihilation oc- 
cured within the bulk BEG after vortices are pulled from 
the chopsticks beams, or if the vortices remained pinned 
to the lasers during the circulation cancellation process. 

To further examine dynamics of circulation cancella¬ 
tion, we performed a numerical simulation of the simul¬ 
taneous return of the beams to the initial position in a 
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return time of = 1.0 s. From the simulation, shown in 
the top row of Fig. 7, we see that the vortices are first 
ripped away from the chopsticks beams within a time 
of tjtr = 0.11. Additional vortices are also seen around 
each beam’s position, but none of these are pinned by the 
beams. Eventually all vortices annihilate one another, 
and the system is left free of any fluid circulation. 

We interpret this process as arising from the laser 
beams now moving against the background fluid flow ini¬ 
tiated by the motion of the chopsticks beams, i.e. the 
fluid flow corresponding to the vortices now pinned to the 
beams. At a nominal time = 0.1 during the beam 
return, the beams are approximately 13.5 /im apart. If 
vortices of opposite circulation are pinned to the beams, 
there is a background flow speed of ^ 52/im/s at each 
beam due to the presence of the vortex that is pinned to 
the opposite beam. If the beams return to their initial po¬ 
sition over distances of 27 fim in 1.0 s, they increase the 
speed of a beam relative to the background flow speed 
by nearly 27 fim/s. The presence of a vortex pinned to a 
beam further increases the relative flow speed at the edge 
of each beam that is closest to the opposite beam (i.e, at 
the boundaries of the channel between the beams). 
Apparently, these conditions are sufficient to rip vortices 
off of the beams, and soon thereafter, to cause the nucle- 
ation of a vortex dipole at each beam, with all vortices 
then annihilating pairwise. 

In a second numerical study of the return of both 
beams to the initial locations, we again performed the 
simulation described above, but first conjugated the BEC 
phase at the end of the chopsticks procedure. Phase 
conjugation is equivalent to reversing the charge of each 
vortex, or to slowly exchanging the positions of the two 
beams by rotating them around their center point, which 
was also examined numerically and is experimentally fea¬ 
sible. After phase conjugation, the beams were returned 
to their initial positions in 1.0 s. Because the beams 
were now moving in roughly the same direction as the 
background flow, the speed of the beams relative to the 
background flow was much less than ^ 52/im/s and 
the vortices remained trapped until the opposite circu¬ 
lations canceled each other when the beams recombined, 
as shown in the second row of Eig. 7. Eurthermore, as 
the third row of Eig. 7 shows, the phase profiles for the 
times at which the chopsticks beams have the equivalent 
spatial separation during the initial forward chopsticks 
sweep are remarkably similar to the return sweep after 
the phase conjugation step. We therefore conclude that 
while the chopsticks procedure cannot be simply reversed 
in order to return the system to the initial state, the pro¬ 
cess does demonstrate reversibility when vortex charges 
are also first reversed. This charge-time symmetry is sim¬ 
ilar to the hysteresis recently reported when spinning up 
or slowing down the persistent superfluid current in a 
toroidally trapped BEC m- 
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FIG. 7. (Color online.) Numerical results of 2D GPE simu¬ 
lations showing the return of both chopsticks beams to their 
initial positions over a return time tr, for beam and BEC 
parameters the same as those of Fig. 3. Phase profiles are 
shown for 61-/im-square regions centered on the midpoint be¬ 
tween the chopsticks beams, with display conditions identical 
to those of the second column of Fig. 3. For ease of com¬ 
parison with the forward sweep simulations shown elsewhere 
in this paper, and specifically for the vortex generation pro¬ 
cedure for a sweep time tg = 0.7 s, shown in the third row 
for sweep conditions identical to those of Fig. 3, the snapshots 
are given from left to right in reverse-time order. The top row 
shows the return of the beams after a 0.5 s hold time, with 
tr = 1.0 s. Beams are moving from right to left, and are in 
the process of merging with one another. At time ijir — 0.07, 
immediately after the beam return starts (top row, rightmost 
image), the vortices are still co-located with each laser beam. 
Soon afterwards, at ijir — 0.11, the vortices are clearly no 
longer pinned (top row, third column). Additional vortices 
are subsequently nucleated (top row, second column) and an¬ 
nihilated, leaving a state with no circulation by ijir — 0.38 
(top row, first column). For the same return sweep, but with a 
conjugation of the quantum phase (second row, also see text), 
the vortices remain pinned to the beams until the beams re¬ 
combine. This return process is nearly identical (aside from 
the phase conjugation) to the initial vortex generation stage, 
shown in the third row for beam separations identical to those 
of the first two rows, but presented for increasing times from 
left to right. 


VI. CONCLUSIONS AND OUTLOOK 


The primary results of this paper are the proof-of- 
principle experimental and numerical demonstrations of 
a deterministic method of creating and manipulating, on- 
demand, two vortices of opposite circulation in a BEC. 
There are numerous parameters involved in determining 
the range of utility of this method, and here we have 


























begun to explore some of those parameters. Our demon¬ 
strations are a step towards the engineering of arbitrary 
vortex distributions in a BEG, and also open up a new 
pathway to precision studies of the interactions of vor¬ 
tices with each other, with sound, with potential barriers 
and obstacles in a BEG, and as a function of other pa¬ 
rameters such as trap geometry and temperature. To 
more fully understand how this technique can be best 
utilized for studies of vortex dynamics and vortex inter¬ 
actions, additional experimental and numerical research 
is required, particularly regarding an assessment of the 
motion of the vortices once the chopsticks beams are 
ramped off, and the ranges of beam powers, sizes, and 
sweep speed that can be successfully used in the chop¬ 
sticks process. 

The goal of creating specified distributions of many 
vortices, on-demand, remains an experimental challenge, 
but one that is within reach. Erom numerical explo¬ 
rations, it appears possible to simultaneously utilize mul¬ 
tiple pairs of chopsticks beams, or to use just two chop¬ 
sticks beams with stationary pinning sites in successive 
vortex creation stages. Additional beam manipulation 
procedures that remove vortices from a BEG can also 
be explored in order to create charge-imbalanced popu¬ 
lations of vortices. These aspects of the chopsticks tech¬ 
nique are numerically explored and presented in a com¬ 
panion article [32] . 
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